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Abstract. We construct an explicit intertwining operator C between the Schrodinger 
group e lt ~ and the geodesic flow g* on certain Hilbert spaces of symbols on the cotangent 
bundle T*Xr of a compact hyperbolic surface Xr = T\D. We also define T-invariant <?*- 
eigendistributions PSj t k,i/j,-u k (Patterson-Sullivan distributions) out of pairs of A-eigenfunctions, 
generalizing the diagonal case j = k in |AZj . The operator L maps PSj^,Vj,-v k to the 
Wigner distribution Wj k studied in quantum chaos. We define Hilbert spaces T-L* PS spanned 
by {PSj^k^j.-Uk}, resp. Hyp spanned by {Wj",.}, and show that £ is a unitary isomorphism 
from Hp S — > The symbols lie in the dual Hilbert spaces. 

1. Introduction 

On a hyperbolic surface Xr = T\D, there is an intimate relation between the spectral 
properties of the laplacian A and those of the geodesic flow g l on the unit tangent bundle 
S*Xr. The Selberg trace formula gives an exact formula for the trace of the Schrodinger 
flow e %t ^ as a sum over closed geodesies, and it may be interpreted as a trace of the pull- 
back operator by g l [G]. Equivalently, eigenvalues of A are (re-parameterizations of) the 
resonances of g l (see [B] for background); see also [Bisl IPo2l |M] for some of the many 
different perspectives on this relation. In this article, we give a yet stronger relation between 
the two flows : we construct an explicit intertwining operator £ (Definition [2]) between 
the Schrodinger flow and the geodesic flow. Our main result, Theorem [HJ is that there 
exist Hilbert spaces of symbols on which £ is a unitary intertwining operator between the 
classical and quantum flow. Much of the problem is to construct the appropriate Hilbert 
spaces, which we denote by ~Hw, 'Hps ( see Definitions 16. HlS~2l) . They cannot be the standard 
Hilbert spaces, L 2 (SX r ) for g*, resp. Hilbert-Schmidt operators for e lt ~, since the spectrum 
of g l is continuous, while that of e lt ^ is discrete; and they are also quite different from the 
Banach spaces constructed in |BT| IBKL] IBL| IFRSt IGL] in the theory of resonances of g l . 

The construction of £, Hw, 'Hps and the proof of the intertwining property grow out of our 
previous work |AZj , where we introduced and studied a family of distributions (that we called 
Patterson-Sullivan distributions) on the unit tangent bundle of a hyperbolic surface. These 
distributions are invariant under the geodesic flow, and we showed that they are closely 
related to the Wigner distributions appearing in the theory of quantum ergodicity. The 
Patterson- Sullivan distributions are naturally constructed from the family of eigenf unctions 
of the laplacian, and we showed that they also arise as residues of dynamical zeta-functions 
at the poles located on the critical line. 
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In this paper, we introduce the family of off-diagonal Patterson-Sullivan distributions, 
and show how they are related to the off-diagonal Wigner distributions (appearing in the 
study of quantum mixing). This construction is a rather straightforward generalization of 
the work done in |AZj . More importantly, we show that these formulae directly lead to an 
operator intertwining the geodesic flow and the Schrodinger flow on the hyperbolic plane 
(or a compact quotient Xr). Roughly speaking, the dual Hilbert spaces W^, resp. T-i* PS , 
intertwined by £ are spanned by the Wigner, resp. Patterson-Sullivan, distributions. The 
main goal of this paper is to construct explicitly this intertwining operator, first on the 
hyperbolic plane, then on a compact quotient, and to investigate some of its properties. The 
existence of this operator opens the way to the construction of a quantization procedure 
satisfying the Egorov theorem in its exact form (without remainder term). 

We have to explain in what sense one can find an intertwining operator between the 
geodesic flow and the Schrodinger flow. The former acts on functions on the (co)tangent 
bundle TXr whereas the latter acts on functions on the base manifold Xp. In fact, we let 
the Schrodinger group act on the space of operators, by conjugation (as in the Heisenberg 
picture of quantum mechanics). Operators have a Schwartz kernel, which is a distribution on 
the product Xr x Xp. By taking the local Fourier transform of the kernel with respect to the 
second component, we get a distribution on the cotangent bundle T*X r , called the symbol 
of the operator. This way, we see that the Schrodinger group acts naturally on the space of 
distributions on T*Xr (in the paper we will always identify the tangent and the cotangent 
bundles by means of the riemannian metric). With this formulation, the Schrodinger flow 
acts on the same space as the geodesic flow, and it is in this sense that we shall intertwine 
their actions. 

1.1. Notation. To state our results, we need to introduce some notation (see §2] for more 
details). We will denote G = PSU(1, 1) ~ PSL(2, R), K = PSO(2, R) a maximal compact 
subgroup, and G/K the corresponding symmetric space, for which we will in general use the 
picture of the hyperbolic disc D = {zGC,|£|<l}, endowed with the riemannian metric 



This is a standard normalization in hyperbolic geometry, but we caution that it differs by 
a constant factor from the normalization used by Helgason [He] ; for us, the L 2 -spectrum of 
the laplacian on D is (— oo, — j]. Hence some discrepancies between some of our formulae 
and Helgason's. 

It is well-known that G can be identified with the unit tangent bundle SJD of the hyperbolic 
disc D (when using the theory of pseudodifferential operators, it is more natural to work on 
the cotangent bundle, but we will always identify both). We will be particularly interested in 
the geodesic flow, which acts on G by right multiplication as follows : for all g G PSL(2, R), 




of the horocycle flow {h u ) u£ -£i, acting by h u (g) = gn u where n, 




G SL(2,R). 
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We will work with two special parameterizations of the unit tangent bundle, identified 
with G. The first one is obtained by writing G ~ (G/K) x K ~ D x B where B is the 
boundary at infinity of D, identified with the unit circle 5* 1 in the Poincare disc model. The 
group K and the boundary at infinity S 1 are identified by the map 

(cos 9 — sin 6 
sin 8 cos 9 

This way, a point in G can be parametrized by the coordinates (z, b), where z G D and 
b G B. In geometric terms, if (z, b) is identified with a unit tangent vector in SD, then b 
represents the (forward) limit point of the geodesic generated by (z,b). The action of G on 
itself by left-multiplication yields an action of G on B (Section EJ). 

We shall also use the following parameterization : denote B^ = {(&', b) G B x B,b ^ b'} 
the set of pairs of distinct points in the boundary. Each oriented geodesic in D is completely 
determined by its (unique) forward limit point b in B and its (unique) backward limit point 
b' 7^ b in B : we will denote 75/ ( & the geodesic going from b' G B to b G B. Thus, B^ can 
be naturally identified with the set of oriented geodesies of D. The elements of G can be 
parametrized by (&', b, r) with (&', b) G B^ and r G R. We identify (&', b, r) with the point 
(z, b) G D x B, where z is on the geodesic 7^, situated r units from the point z^y G 7^ 

closest to the origin o = f eK G D. 

Our final goal is to obtain formulae that are valid on a compact quotient of D; that is, we 
consider a co-compact discrete subgroup T C G. We assume it has no torsional and contains 
only hyperbolic elements. Then the quotient Xr = T\D is a compact hyperbolic surface. 

1.2. Quantization, Wigner distributions and Patterson-Sullivan distributions. A 

quantization procedure adapted to the hyperbolic disc was defined in [Z3j . using Helgason's 
version of the Fourier transform |Hej . For (z,b) G D x B, define the Busemann function 
(z, b) as the signed distance to o of the horocycle going through the points z G D, b G B. The 
family of functions z 1— > e^ +ir '^ z,b ' (r > 0,b G B) forms a basis of generalized eigenfunctions 
of the laplacian on L 2 (D) [HeJ. The hyperbolic pseudodifferential operators introduced by 
|Z3] are defined by 

(1.1) Op(a)e^ +ir){, ' b) = a(;b,r)e^ +ir) ^ b \ 

if a = a(z, b, r) is a function on D x B x R ~ S*D x R which must have "reasonable" decay 
and smoothness properties (section Hj). The function a is called the symbol of the operator. 
We note that by choosing r < instead of r > we obtain another basis of generalized 
eigenfunctions of the laplacian. We also note that the Schwartz-kernel of the operator is 
formally given by 

(1.2) K a (z,w)= [ [ a(^,6 / ,r)e ( 5+^< 2 ^e ( 5- ^r )< ^t '^dp(r)d6 , , 

where dp(r) is the Plancherel measure defined in §3.11 Paley- Wiener type theorems relating 
the decay and regularity of a and those of K a will be recalled in §4.21 In most formulae we 




This assumption is probably not necessary. 
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must assume that K a decays sufficiently fast away from the diagonal {z = w}. This implies 
in particular that the corresponding symbol a has a holomorphic extension to r G C. 

Let Xr be a compact quotient of D as above, and fix an orthonormal basis (<f) k ) of L 2 (Xp) 
formed of eigenfunctions of the laplacian. We use the standard notations in hyperbolic 
spectral theory : the eigenfunctions 4> k satisfy 



where v k = ir k G R U ?R is called the spectral parameter (on a compact surface, only a 
finite number of r k s are imaginary). For each eigenvalue there are two possible choices for 
the spectral parameter. 

If a is T-invariant, it was shown in \Z3\ [Z] that Op (a) preserves the space of T- invariant 
functions. We denote by Op r (a) the operator Op(a) acting on T-invariant functions. The 
'Wigner distributions^" Wj^ are defined on ST) x R ~ G x R by the formula 

(1.3) / adW jtk = (Op(a)0 i ,0 fe ) D , 

iSDxR 

for a a function on SY) x R, with appropriate growth and smoothness properties. The 
distribution Wjk is invariant by the action of T on SY), and thus can be used to define a 
distribution Wj k on the quotient SXr x R ~ Y\G x R : if a is a function on Y\G x R, in 
other words a T-invariant function on G x R, we define 



/ a dWj k = xa dW jtk , 

Jsx r x~R Jsdxh 



where \ is a smooth fundamental cut-off function for the action of Y (see §3.3p . It can easily 
be seen that this definition does not depend on the choice of %. Moreover we have 



;i.4) / adWj ]k = (O Pr (a)^-,0 fe ) L 2 (Xr) = Tr Xr Op r ( 



a) o {(pjd$(p k ), 



for a a function on S"Xr x R, with appropriate growth and smoothness properties (see Section 
H]for a detailed discussion). 

The Wigner distribution Wj k may also be expressed in terms of the boundary values of 
the eigenfunctions <f>j, <f> k . The boundary values T k)Vk of (p k is a distribution on the boundary 
B, with the property that 

M*) = [ e^^T k>Vk (db), 
Jb 

for all z G D. It depends on the choice of a spectral parameter u k , and is unique if we pick 
u k such that | + v k ^ 0, -1, -2, • • • ([He, Theorems 4.3 and 4.29]; see also [He2j ). With a 
slight abuse of notation we shall denote T Vk instead of T k ^ Vk . Using the definition of Op, we 
have 

/ adWl k = [ X (z,b)a(z,b,r j )^e^ + ^ b %(db)Vol(dz) 

J5X r xR JUxB 

where \ is a smooth fundamental cut-off, see §3.31 



2 This terminology is normally used in a euclidean context, but it is sometimes convenient to extend it to 
non-euclidean geometries. 
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It follows from its definition that Wj k is an eigendistribution of the quantum evolution. 
Define 

(1.5) a'(Op(a)) = e" ft T Op(a)e ft T. 
We then have 

(1.6) (a\0 Vv {a))fa,fa) = e^-^^iOpM^Ak) = (Op r (a)fa, fa) 

where the last identity holds only for real values of rj, r k . We henceforth denote by V 1 the 
operator on symbols, defined formally by 

(1.7) a*(Op(a)) = Op(y*(a)). 

See Section H] for precise assumptions on a. We shall denote Vp when we want to stress the 
fact that V 1 act on T-invariant symbols. 

The Wigner distributions have been studied a lot in the theory of quantum ergodicity and 
quantum mixing. In the context of T\C7, the Wigner distributions of this article are studied in 
[Z| [Z3| IW] ISVj IEV2J as well as in |AZj . In this paper we introduce the family of (off-diagonal) 
Patterson- Sullivan distributions. They are also constructed from the eigenfunctions <pj. Take 
the boundary values T^ v .{db) of (f)j and T k _ Uk (db) of fa. 

Definition 1. PSj : k,uj-u k (db' , db, dr) is the T -invariant distribution on B^xH ~ G defined 

by 

pc (Ail Ah a \ T u 3 {db)T_ Uk (db') , +W)T 

PS {j , UjUk ,- Uk) (db,db,dT) = ^_ ._- e ( > + ^ dr. 

We note that the Patterson-Sullivan distributions depend on the eigenfunctions fa, fa, but 
also on the choice of the spectral parameters Vj, v k (in contrast with the Wigner distributions, 
which depend only on the eigenfunctions); hence the notation PS(j, V j),(k,-v k )- 111 the sequel, 
we will in general use the shorter notation PS Vj) - Vk , although it is slightly abusive. In 
Proposition 16. 1[ we check that the distributions PS v .- Vk are (right)-P-invariant distributions 
on x R ~ G. Besides, since the geodesic flow reads 

( ? t (6',6,r) = (6 / ,6,r + t), 

they are eigendistributions for the geodesic flow in the sense that 

(1.8) 9\PS Vj ,- Vk = e-^+^PS^ = e a ^-^PS^_ Uk 

(the last identity holds only for real values of rj,r k ). As a result, PS v -- Vk induces an 
eigendistribution PS'E. _ v of the geodesic flow on r\C7 = SX. r , defined by 

(1.9) / adPS T = [ { X a)dPS Vj ,- Vk , 

JT\G JG 

for every smooth T-invariant function a. Once again x is a smooth fundamental domain 
cutoff, see §3.31 When uj + Vk = 0, these Patterson-Sullivan distributions are the family 
of diagonal Patterson-Sullivan distributions, introduced in |AZ] (these are invariant by the 
geodesic flow). Recently, Hilgert and Schroder have extended the definition and properties 
of off-diagonal Patterson-Sullivan distributions to more general symmetric spaces |HilSc[ 
ISchDiss] . 
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It was pointed out in [Z] that the distribution e u . {z, b)e^ Yo\(dz)db d = e^ +u ^ b) T Uj {db) Vol(ds) 
on SXr is a joint eigendistribution of the horocycle and geodesic flows, contained in the 
(spherical) irreducible representation of G generated by <pj. The PS-distributions are new 
objects, and it is illuminating to express them in terms of these more familiar ones. 

Proposition 1.1. The Patterson-Sullivan distributions are given by the (well-defined) prod- 
ucts [e Vj ■ t-e-u k ) ■ 

PS v .- Vk {dV, db, dr) = — (e v . ■ LtZ^) (z, b) e {z ' b) Vo\(dz)db, 

Z7T 

with (z,b) ~ (b',b,r) and where i is the involution (V ,b,r) \- > (b,b',—r), corresponding to 
(x,£) i — y (x, — £) on SD or to the action of the non-trivial element of the Weyl group, g i— y gw 
-1 

on G, w 



1 



Remark 1.1. We note that the Wigner distributions are naturally defined on S'Xp x R, 
whereas the Patterson-Sullivan distributions were only defined on SXr- In order to relate 
both, we need to extend the latter to SXp x R. We do so by identifying SX r with SX r x 
{ Vj 1 2 rk }. In other words, we extend the PS'-distributions to SX r x R by considering 

2 

at least for real values of r J; r^. Doing so, we must pay special attention to the case of low 
laplacian eigenvalues, when r is imaginary. For imaginary rj,rf., the formular above will be 
generalized to 

PSl QS^-vt, 

J 2i 

and our formula will hold for functions a(z, b, r) that extend holomorphically to imaginary 
values of r (corresponding to kernels decaying fast enough away from the diagonal). 

1.3. Definition of the intertwining operator. 

Definition 2. The intertwining operator 

C : C C °°(G x R) -> C(G x R) 

is defined by 

(1.10) Ca(g, R) = J (1 + u 2 )~^ +iR) a{ga T log( i + „2 ) w u , r) e - 2i(R - r)l 'drdudr. 

Here, n u G N is the one-parameter unipotent subgroup whose right-orbits define the horo- 
cycle flow and a t 6 A is the one parameter subgroup whose right-orbits define the geodesic 
flow (see AO) . 

Extend the geodesic flow to G x R by the formula 

(1.11) G t (g,r) = (ga rt} r). 

We will also consider the geodesic flow as an operator acting on functions, by composition : 
for a function a on G, we denote g t a = aoj f , and for a function on G x R, G l a = a o G 4 . 
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1.4. Statement of results. The main result of this paper is that we have the intertwining 
relation 

(1.12) CoV t = G t o£. 

We prove f ll . 12[) in several levels. At first, we prove that the dual (or adjoint) equation holds 
if we apply both sides to the Patterson-Sullivan distributions. At this level, we do not specify 
the domain and range of the dual intertwining operator £ on symbols. Indeed, the domain 
and mapping properties of C on Xp are rather subtle, and we only introduce the domain in 
§0 after first studying the mapping properties of C in §5] on the hyperbolic disc (Theorem 
[7]) . The final result is Theorem |HJ 

We now state the first level of the result, whereby the intertwining relation comes from 
an exact relation between the Wigner and Patterson-Sullivan families of distributions. If a 
is a function on G, decaying fast enough, and v G C, define 

L u a{g)= i (l + u 2 y^ +u) a(gn u )du. 
Jn 

If a is a function on G x C, and v e C, define the function a u on G by a u (g) = a(g,u). 
The following (somewhat imprecise) statement asserts that the adjoint of £ maps Patterson- 
Sullivan distributions to Wigner distributions. It is stated precisely in §6.31 

Theorem 3. Let a = a(z,b,r) be a T -invariant function, with 

a(z,b,r) = J^a(7 • 2,7 " & > r )> 
7er 

with a satisfying adequate decay and smoothness properties. Then we have 
Wj^a) = 2 1+ ^PS u ^ Uk {L^~a U] ) = PS^ k (Ca). 

Theorem [3] is first proved on the hyperbolic disc (Section [5]) by introducing analogues 
of the Wigner and Patterson-Sullivan families there, and then goes down to the quotient 
Xr (Section [6]). The "adequate decay and smoothness properties" are described in §6.31 

We now state the final result in a somewhat imprecise form. In §6.3[ we define Hilbert 
spaces and Wp S . The Wigner distributions form an orthonormal of Ti^, while the 
Patterson- Sullivan distributions form an orthonormal basis of T-L* P g- Both are modelled on 
the Hilbert space of Hilbert-Schmidt pseudo-differential operators. The dual Hilbert spaces 
Tiw, resp. Hps are the spaces of symbols on which C acts as an intertwining operator. The 
space HS^ 3 is defined in Definition 14.11 

Theorem 4. The interwining operator '■ H*ps ~ >* H-w ^ s an isometric isomorphism, and 
£rj sends PS Vj ~ Vk to Wj^- Dually, there exists a class of automorphic symbols, a 6 , 
so that 

C r o V^a = G r o C r a, 
as an equality between two elements of Hps- 

There are similar partial results for the wave flows e' lty ^^ and the modified wave flow 
eV- A -i (see ED - 
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1.5. Asymptotic equivalence of Wigner and Patterson-Sullivan distributions. In 

[AZJ, it is proved that, after suitable normalization, the diagonal Wigner distributions and 
Patterson-Sullivan distributions are asymptotically the same in the semi-classical limit. The 
same is true for the off-diagonal elements: 

Theorem 5. Let a e C°°{T\G). Given a sequence of pairs (vj n ,Vk n ) of spectral parameters 
with —ivj n — > +00 and \vj n — Vk n \ < To f or some r > ; we have the asymptotic formula 

[ a(g)Wl kn (dg) = 2^^ f^-Y' \-S [ a(g)PS^_ u Jdg) + 0{y£). 

JSXr \ r k n J Jsx r 

The proof is very similar to that in the diagonal case in [AZJ, starting from Theorem [3] 
Hence we only sketch the key points in §HJ This result has been extended to more general 
symmetric spaces by Hilgert and Schroder (HilScJ. 

In |Z2] (see also |Z5j). it is shown that the off-diagonal Wigner distributions Wj n ^„ with 
jn 7^ k n and with a limiting spectral gap r Jn — r kn — > r tend to zero when the geodesic 
flow is mixing, at least after the removal of a subsequence of spectral density zero. It then 
follows from Theorem [5] that: 

Corollary 6. Take a sequence of pairs (j n , k n ), with j n 7^ k n and r Jn — r kn — > r . Assume 
that this sequence has positive density, in the sense that 



lim inf 



UK 1 




|<A} 


HO", 


Vi\ 


<A} 



> 0. 



Then there exists a subsequence of full density such that r k J PS Vn ^ Vk — > 0. 

1.6. Relations to other work. The existence of the intertwining operator is rather unex- 
pected from the viewpoint of microlocal analysis and quantum chaos, but is quite natural 
from the viewpoint of automorphic distributions and invariant triple products [BR] [BR2, 
IMSj iDf ISVj ISV2] . where it may be interpreted as intertwining the family of Wigner triple 
products £ w (a,4>j,(f> k ) = (Op(a)0j, <f) k ) and the family of Patterson- Sullivan triple products 
£ PS (a, <pj, 4>k) = (a,dPS Uj - Uk ). It follows from general principles that there exist constants 
C r . trk £ PS (a,<j)j,(f) k ) = C rjtrk £ w (a,(j)j,(f) k ) and essentially £ is an integral operator with ma- 
trix elements C rjj1 . k . Explicit formula are given in |AZ2j relating the triple products when 
evaluated on automorphic symbols of a fixed weight. This approach through explicit formu- 
lae and representation theory is a (less global) alternative to the study of C. But C might 
have an independent interest in representation theory as the intertwining operator between 
these two families of triple products. 

The relations between Wigner and Patterson- Sullivan distributions, and the exact formulae 
relating them in |AZ2j . shed some light on the limit formula for quantum variances of Wigner 
distributions proved by Luo-Sarnak |LSj and Zhao [Zh] . The quantum variance for a zeroth 
order pseudo-differential operator A is defined as 



(1.13) 




(A(A) = #{j : A, < A}) 
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where J a du is the Liouville average of the principal symbol a of A. It was suggested by 
Feingold and Peres [FPJ that the quantum variance should tend to the following way : 

The bilinear form B should be proportional to Pa o ,a o (0), where 



(1.14) P4>,ip(t) = / (p(x)4'(g t x)duj(x) — / 4>duj / tpduj 

is the "dynamical correlation function", and p^^ij) = e~ %Tt p^(t)dt is its Fourier trans- 
form. In |LS[ IZh] a version of this conjecture (with additional arithmetic factors) was proved 
for the basis of Hecke eigenfunctions. 

We postpone further discussion to |AZ2j . but would like to state a conjecture that con- 
nects Patterson-Sullivan distributions to the quantum variances in |LSt IZh] . Unlike Wigner 
distributions, Patterson-Sullivan distributions are defined independently on the choice of a 
quantization procedure. In view of Theorem^, it is natural to consider the classical variances 
for the diagonal PS-distributions: 



[1.15) PS a (X) d ^ ^ K^SW - J adu\ 



Here, PS Vj - Vj = f — PS Vj - Vj are normalized PS-distributions (see |AZ] ) so that 

( 1 , PSvj , - uj ) sx r 

the statement is correct for constant functions. Since they are sums of g*-invariant bilinear 
forms, they should have a closer relation to dynamical correlation functions than variances 
for Wigner distributions. 

Conjecture 1. Let Xr be a compact hyperbolic surface. There exists a constant Cx r such 
that 



PS a {\) 



Cx r Pa o ,a o (0) 



Vx ' 

We note that the quantum variance associated with Patterson- Sullivan distributions may 
still depend on the choice of a basis of eigenfunctions. 

Finally, we point out a possibly tenuous relation of our intertwining problem to the one 
studied by Bismut on locally symmetric spaces of non-compact type in Chapter 10 of [Bis]. 
On the infinitesimal level, we are intertwining the generator of the geodesic flow to the 
operator P taking a symbol a to the symbol of [Op(a), A]. As discussed in [Z3] the latter 
operator P has the form H 2 + 4X^ + irH, where H 2 + is elliptic along the stable 
foliation and irH is the semi-classical operator of order 2 where H generates the geodesic 
flow. By comparision, Bismut 's hypoelliptic laplacian is essentially the weighted sum 
of the harmonic oscillator on the fiber of T*Xp and bH. We note that the H terms are 
identical if we set b = r, while the other terms are in a sense orthogonal (Bismut's is vertical 
while ours is horizontal). But both have the essential property that as the semi-classical 
parameter b = r — > oo, the operators converge to the generator of the geodesic flow. There 
is a possible parallel of our conjugation problem to the conjugation between the hypoelliptic 



laplacian and a certain elliptic operator in [Bis] (Chapter 10). We encounter similar problems 
in defining the domain of C and its inverse. 
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Notational issue : In what follows we have to face the issue that we are sometimes 
using the sesquilinear pairing between two L 2 functions (or more generally, two elements of 
a complex Hilbert space), and sometimes the bilinear pairing between a distribution and a 
test function (more generally, an element of a vector space and a linear form). We will try 
to keep distinct notations to avoid confusion, denoting by (a, b) the scalar product of two L 2 
functions (linear w.r.t. a, antilinear w.r.t. b), and by T(a) the pairing between a distribution 
T and a test function a. More generally, we shall use the bracket notation (., .) exclusively 
for sesquilinear pairings. 

If L is a linear operator on a space endowed with a sesquilinear form, we shall denote by 
U its adjoint in the hermitian sense, that is, (a, Lb) = (L*a,b). 

We will denote L$ the adjoint of L in the usual sense of duality : if T is a linear form, 

then L<fT is the linear form defined by (LfT)(a) = T(La). Thus, L i— > is linear whereas 
L i— > U is antilinear. 

If T is a distribution and <3> a diffeomorphism, we shall also denote by $j)T the pushforward 
of T by $ : ($j|T)(a) == T(a o $). The two notations should not interfere. 

2. Coordinates on SJD. 

2.1. Dynamics and group theory of G = PSL(2,~EL). A set of generators of the Lie 
algebra si (2, R) is given by 

/ 1 \ / 1 \ / -1 

h=[ , x + =l , y = 

\ -1 / \ / \ 1 

The subgroups they generate are denoted by A, N, K respectively. We also put X_ = 

\ _ 
, and denote the associated subgroup by N. In the identification SD = PSL(2, R) 

1 Oj 

the geodesic flow (<7*)teR is given by the right action of the group A of diagonal matrices 

/ e*/ 2 \ 

with positive entries : g H- ga t where a t = I . The action of the horocycle 

V e"*/ 2 / 

flow (/i u ) nS R is defined by the right action of N, in other words by g i— > where n u = 
1 u\ / 1 

. We shall also denote n u = I 
1/ \ u 1 

2.2. Adapted coordinates. As explained in the introduction, the identification G ~ 
(G/K) x (G/NA) ~ (G/K) x X leads to the coordinates (z, 6) (where « G D, 6 6 G/NA ~ 

~ S* 1 ) to parametrize points in G ~ ST). The identification G ~ (G/K) x (G/NA) is G- 
equivariant, and thus the action of g on itself by left-multiplication reads g-(z,b) = (g-z,g-b), 
where on the first component G acts by isometry on the symmetric space G/K, and on the 
second coordinate G acts on the boundary G/NA ~ K. 

We denote by o the origin eK in G/K, and by (-2,6) the signed distance to o of the 
horocycle through the points z E~D,b E B. This notation follows [He] , but we warn again 
that our normalization of the metric differs by a factor 2 from Helgason's. 
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We also use the identification G = ST) = B^ x R. It is based on the identification of 
£?( 2 ) = B x B\A (where A C B x B is the diagonal) with the space of oriented geodesies 
of D. To (b', b) G B x B \ A there corresponds a geodesic 75/ & whose forward endpoint at 
infinity equals 6 and whose backward endpoint equals V . The choice of time parameter is 
defined so that (b',b, 0) is the closest point zy )6 to the origin o on 7 6 / j6 , and (b',b,t) denotes 
the point t units from (&', 6, 0) in signed distance towards b. We note that (z&y, 6) = ( &')• 
We define g>(&', b) G PSU(1, 1) to be the unique element satisfying 
<?(&',&)• 1 = 6, 

• (-1) = fo', 

• o = z v , b , 

where 1, —1, b, b' are points of the boundary B = S 1 , seen as the unit circle in C in the disc 
model. We thus identify B^ x R ~ G by (b',b,t) h-> g(b',b)a t . In these coordinates, the 
action of g G C7 by left-multiplication is expressed by 

/,/ , ,x / w h , . (9 ■ o,g -b) - (g ■ o,g -b')\ 

(2.1) g- {b ,b,t) — \g-b,g-b,t+ J. 

We will need the following formula : 
Lemma 2.1. 

|o-6'| 

log — h (5(6 , 6)a r ^ M • o, 6) = r, 

where b,b' are seen as elements of S 1 C C ; and |6 — &'| t/iezr usuai distance in C. 
Proof. To prove this, we use the identity 

(ga t n u ■ o, g - 1} = (a t n u • o, 1) + (# • o, g - 1) = t + (g ■ o, g - 1) 
to reduce the lemma to the claim that 

(2.2) ( ^ )6 ). 0) ^ )6 ).i) = _ log t^l. 
However, a basic identity gives 

(2.3) \g0 - g p>\2 e (9-o,g-P)+(g-o,g-P') = \q _ £>|2_ 

If we let (3 = l,/3' = -1 (so that 5(6', 6)1 = b,g(V,b)(-l) = b') and recall that (g(b' ,b)-o,b) = 
(g{b',b) ■ o,6'), then (E3D implies 

(2.4) 4 = \b - b'\ 2 e^ Zb '' b,b)+{ - Zb - b '' b '^ = \b- b'\ 2 e 2{ - Zb ' - b ' b \ 
which completes the proof of (12.21) and hence of the lemma. 

□ 

2.3. Time reversal. Time reversal is the map 1 : (x,£) — > (x, — £) on the tangent bundle. 
In the coordinates (b',b,t) it takes the form, 

(2.5) .(&',M) = (6,6',-t). 

That is, it reverses the endpoints of the oriented geodesic 7^ and preserves the point z^y 
closest to o. In the group theoretic picture, time reversal is given by the action of the 

non-trivial element of the Weyl group, g 1— > giv on C7, where w = 
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2.4. A coordinate change. The formulae below are useful at several places in the paper. 

(2.6) n u = A: u a-iog(i+u 2 )"-/(«), 

where f(u) = and where 

This comes from the explicit calculation, 

J u _\ / 1 o 

(2.7) 










: 








This formula implies that the element g = n u a t G G corresponds to the endpoints 6 = 1, 
b' = k u w G K, in other words, using the S* 1 model, 

b' = e 2W , where e td = , U + 

1 /-I o 



By calculation, we find |6' — 6| 2 = and db' = These calculations also show that 

(n u a t , 1) = t 

and 

(n u a t ,b'} = -t + log(l + « 2 ). 
For t = log ( 1+ " \ we see that (n u a t ,l) = (n u a t ,b'), and thus n u a losil+u 2 } = g Vh (with 6 = 

2 

1,6' = e 2l6> as above). It follows that g = n u a t G G has the coordinates (6', 6, r) = (e 2i6> , 1, t — 

lQg(l+M 2 ) \ 

2 /' 

3. Harmonic analysis on the hyperbolic disc and its compact quotients 

3.1. Poisson 1-form, Haar measure and Plancherel measure. We shall denote by 
db the normalized Haar measure on K, identified with the boundary B or with S 1 . The 
Poisson 1-form is defined by 

(3.1) P(z,b)db = e {z ' b) db. 
Using the identities 

(3.2) (g-z,g-b) = (z,b) + (g-o,g-b), 
and 

(3.3) ±g - b = , 
it follows that 

(3.4) P(g-z,g-b)d(g-b)=P(z,b)db. 

Haar measure on G is denoted dg. In terms of z, 6 coordinates it is given by 

(3.5) dg = P(z,b)Vol(dz)db, 



GEODESIC FLOW AND SCHRODINGER FLOW ON HYPERBOLIC SURFACES 13 

where Vol(dz) is the hyperbolic area form. Under the identification G ~ SD, the Haar 
measure on G is the same as Liouville measure on ST). In the (b',b,t) coordinates, Haar 
measure reads as follows: 

Lemma 3.1. Under the identifications (b,b',t) — g — g(b,b')at — (z,b), we have 

db 6?) db' 

dg = P{z, b) Vol{dz)db = 4tt " _J^ 2 ® dt. 

The fact that the measure njzmz ® dt is invariant under the action of G follows from 
formulae (12. ip . (I2.3p . (13. 3p . We leave aside the calculation of the normalization factor, which 
can be done thanks to the formulae in §2.4[ but does not play a very important role. 

Non-Euclidean Fourier analysis is based on the family of non-Euclidean "plane waves" 

e^{z)^e^^ b \ 
v G C, b G B. They are complex-valued eigenf unctions of the laplacian : 

Ae„ i6 = - Q - v 2 ^ e v>b . 

The L 2 spectral decomposition of the laplacian on D only requires the tempered spectrum, 
that is, the case v = ir where r G R (corresponding to a laplacian eigenvalue \ — v 2 > |). 
The Helgason-Fourier transform of a function / on D is defined by 



Ff(b,r)= / e^- ir ^f(z)Vol(dz) 



D 

b G B, r G R. The Fourier transform automatically has the following symmetry property 
under r i— > —r : 

(3.6) I Ff(b,r)e^ +ir){z > b) db= [ Ff(b,-r)e^- ir){z > b) db, 
Jb Jb 

for all z G D and r G R. Plancherel measure is the measure on R defined by 

(3.7) dp(r) = — r tanh(7rr)c?r, 
and the Fourier inversion formula reads 

f{z) = - I I Ff(b,r)e^ +ir){z > b) dp(r)db= I I J -f(b,r)e^ +ir){z > b] >dp(r)db, 

^ J R, J B J R,_j_ J B 

see [He]. We have the Plancherel formula for / G L 2 (D), ||/||l2(d,Vo1) = ll-^/IU 2 (BxR +l( i6xdp(r))- 

3.2. Integral representation of eigenfunctions. We now consider Fourier analysis on 
the quotient X r of D by a discrete co-compact subgroup T C G. 

Theorem 3.2. ( [Hej . Theorems 4-3 and see also |He2j ) Let be an eigenfunction with 
exponential growth, for the eigenvalue X = — Q — z/ 2 ) G C. Then there exists a distribution 
T v ^ G T>'(B) such that 

Jb 

for all z G D. The distribution is unique if | + v ^ 0, — 1, —2, • • • . 
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The distribution T v $ is usually called the boundary values of cf> (for the spectral parameter 
is), in analogy with the theory of boundary values of harmonic functions. This theorem 
applies, in particular, to a r-invariant eigenfunction of the laplacian, since such a function 
is bounded. By uniqueness of T u ^, we see that being r-invariant is equivalent to 

(3.8) l^Tu^db) = e-^ + ^°^T u ^db) 
for 7 G T and b G B. 

3.3. Fundamental domains and cutoffs for T\G. We denote by D a fundamental 
domain for the action of V on D = G/K. We use the same notation for the fundamental 
domain D lifted to G. 

When dealing with integrals against irregular distributions, it is convenient to replace the 
characteristic function of a fundamental domain by a smooth (compactly supported) cutoff 
X on G satisfying 

n* = i, 

where we define the periodization operator IT by 

n x(#) = ^2x(ig)- 

■yer 

Existence of such functions x is obvious. We will call such a function x a smooth fundamental 
cutoff for the action of T on G. When needed, we may assume that x is a right-fT-invariant 
functions, that is, x{ z ->b) = x( z )- 

Let X-, x' De two smooth fundamental cutoffs. We will use repeatedly the following : if T 
is a r-invariant distribution on G, then T(xf) = T(x'f), for any / G C°°(T\G) (seen as a 
T- invariant function on G). To see this, write 

(3.9) T( X f) = T( X f.(U X ')) = T(U( X f). X ') = T(fx'). 

4. Pseudo-differential calculus on the Poincare disc 

Throughout this article, we use a special hyperbolic calculus of pseudodifferential operators 
introduced in [Z3j . In the hyperbolic calculus, a complete symbol a(z,b,r) ((z,b) G DxB,r G 
R) is quantized by the operator Op (a) on D defined by 

Op(a)e ir<b (z) = a(z,b,r)e^ +ir ^ 

for z G D, b G B and r G R+. By the non-Euclidean Fourier inversion formula, we define 
Op(o) on C C °°(D): 

Op(a)u(z) = [ [ a(2,6,r)e ( 5+ ir )^^u(6,r)dp(r)d6. 
Jb Jn+ 

We recall that the measure dp(r) = tanh(7rr)<ir is the Plancherel measure for G ( 13. 71) . 
At the formal level, the kernel of Op (a) is thus given by 

(4.1) K a {z,w)= I [ a(z,b,r)e^ +ir){z ' b) e^- ir){w ' b) dp(r)db. 

Jb Jri 
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Now assume that a has the following symmetry w.r.t. the transformation r i-» — r: 
(4.2) / a(z,b,r)e^ 2+tr ^e^ 2 - ir ^db= [ a(z, b, -r)e^ 2 ~ ir ^ e^ 2+ir ^ b > db 



for all z, w G D and r G R. It then follows from the Plancherel formula for the hyperbolic 
Fourier transform that we can recover the symbol from the kernel by 

(4.3) a(2,6,r) =e~<5+ < r><*. 6 > I K a (z, w )e^ +ir){w > b) Vol(dtw) 

Jd 

for all r G R. In this case, formula (14. ip holds with J" R replaced by J R . We now discuss 
several particular classes of symbols a. 

4.1. Hilbert-Schmidt operators on D and L 2 symbols. We recall that if H is a Hilbert 
space, the algebra of Hilbert-Schmidt operators on % is the algebra of operators A for which 

the trace Tr is finite; it is endowed with the inner product (A, B) m = Tr ABl It is well 
known that the Hilbert-Schmidt operators on H = L 2 (M, dv) for any measure space form a 
Hilbert space isomorphic to L 2 (M x M, dv x dv). In the case M = D and v = Vol, we will 
denote HS(D) the space of Hilbert-Schmidt operators on D. We denote L^(C7xR, dgxdp(r)) 
the space of functions in L 2 (G x R, rig x dp(r)) that have the symmetry (14. 2 p with respect 
to the Weyl group. We endow it with the norm 



\a(z, b, r)\ 2 P{z, b) Vo\(dz)dbdp(r). 



\a\\\ 2 = - I [ [ \a{z,b,r)\ 2 P{z,b)Vo\{dz)dbdp{r) 

w 2 y D J B J R 

J D JB JR + 

The following is a consequence of the Plancherel formula : 
Proposition 4.1. The quantization map a — > Op(a) defines a unitary equivalence 

L 2 W (G x R, dg x dp) ~ HS(D). 

In other words 

12 _ 1 / / / I, ., _ 7. ,.M2 



l|Op(a)||^ s(D) = - / / \a(z } b,r)\ 2 P(z,b)Vo\(dz)dbdp(r). 

Z JT> JB JR 

It is then clear that the time evolution V 1 (11.71) defines a unitary operator on L 2 V {G x 
R, dg x dp) . 

4.2. Schwartz class and associated symbols. Schwartz functions on G were first defined 
by Harish-Chandra [HC66J; the definition was extended to G/K by Eguchi and his collabo- 
rators |Eg74 , Eg79|. Writing the hyperbolic disc as G/K, f belongs to the Schwartz space 
C P (G/K) (for < p < 2) if and only if / is a function on G which is right-i^-invariant, and 

svL-pVo(gK)- 2 ^(l + d(gK,o)y\LRf(g)\ < +oo, 

for any q > 0, and for any differential operators L, R on G which are respectively left- and 
right-invariant. Here (p Q stands for the spherical function on G/K, (p (z) = f e^ z,h 'dh. It 
satisfies tp (z) x d(z, o)e~ d( - z '°^ 2 as the hyperbolic distance d(z, o) — > +oo. Functions on 
C P (G/K) are, in particular, in LP (they are sometimes called Schwartz functions of L p -type). 
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The Fourier transforms of Schwartz functions of L p -type were characterized by Eguchi 
Eg74l |Eg79] : letting e = e(p) = f - 1, F (C P (G / K)) coincides with the space C(B x R e ) w 



of functions u on B x R such that 

- u extends holomorphically to the strip R e = {|Qm(r)| < §} (this condition is empty for 
6 = 0); 

- on this strip (or, in the case e = 0, on the real axis), we have a bound 



(4.4) sup(l + \r\) q 

(b,r) 



Du(b, r) 



dr a 



< +oo, 



for all q > 0, every integer a, and every K- left-invariant differential operator D acting on B 
(here we use the identification B ~ K); 

- besides, u must satisfy the symmetry (13. 6 p (this symmetry condition with respect to the 
Weyl group is indicated by the subscript w)- 

It is clear from this characterization that the Schwartz space C P (G/K) is stable under 



We now define the space JC p ' g (G/K x G/K) (resp. K VA (G/K x G/K), !C P (G/K x G/K), 
fC q {G /K x G/K)) of kernels of operators sending C P {G/ K) continuously to C q {G /K) (resp. 
(C P (G/K))' to {C q {G/K))', C p {G/K) to {C q {G/K))', {C P {G/K))' to C q {G/K)). We denote 
the corresponding symbol classes by S p,q (G/K x BxR)^ S Piq (G/K x B xR) w , S P (G/K x 
B x H)w and so on. 

All these classes of operators are obviously stable under conjugation by e lt ~ , and thus V 1 
preserves the corresponding symbol classes. 

We will in particular consider the space K™{G/K x G/K) of "smoothing" operators, 
sending (f]C p (G/K))' to f] C P (G/K). Corresponding symbols a(z, b, r) are characterized 
by the fact that a{z,b,r)e^ +tr ^ belongs to P\ e P\ p C(B x R e ;C p (G/K)) w (i.e. functions 
a(z,b,r) with the C(B x R e )-regularity in the (b,r) variables, taking values in C P (G/K)). 
We will denote this space of "smoothing" symbols by 

(4.5) := S™(G/K xBxR) w := flfl^ X ^C P {G/K)) W . 

e p 

4.3. Op(a) and Op r (a). A key point of the non- Euclidean pseudo-differential algebra is that 
it is automatically left invariant |Z3j . We say that a symbol a is T- invariant if a(j-z, 7-6, r) = 
a(z,b,r). Denote t g the action of g G G on functions on G/K, defined by t g f(z) = f(g~ l z). 
We recall from [Z3J that a being T-invariant is equivalent to having [t y , Op(a)] = for all 
7 G T. This commutation relation is also equivalent to the fact that K a (j-z, j-w) = K a (z, w). 
In this case, one may view Op(a) in either of two ways: as an operator on C C (D) or as an 
operator on T-invariant functions. The operators differ in the domains they are given. When 
we need to emphasize that the action takes place on T-invariant functions we denote the 
operator by Op r (a). More generally, if an operator A defined on functions on D commutes 
with all the translations £ 7 (7 G T), we shall denote by A r the same operator acting on 
T-invariant functions. 

By the decay properties of the spherical function ip , we see that L 2 (Xp) can be continu- 
ously embedded in (C P (G/K))' if p < 1. As a result, if the kernel K(z, w) is T-invariant, and 
is such that \{z)K{z, w) G K P {G/K x G/K) (where x is our smooth fundamental cut-off 
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§3.3p . then K defines naturally a bounded operator on the quotient, L 2 (Xr) — > L 2 (Xr) : 
for G L 2 (Xr), one can define Kip by the identity 

(K<p,tfj)x r = (xK<f),ifj) D , 

for all if) G L 2 (Xr) . Besides, this definition does not depend on the choice of the fundamental 
cut-off x- 

We can rephrase this in terms of symbols. Assume that a(z, b, r) is a T-invariant symbol, 
and is of the form a(z,b,r) = X] 7 er ^(7 ' z il ' ^i r ) where a G S£ and p < 1 (this holds, 

for instance, if a(z,b,r) = f x(z)a(z, b, r) belongs to Sg). Then we can define the bounded 
operator Op r (a) : L 2 (X r ) — > L 2 (X r ) by 

(Op r (a)0,^) Xr d = (Op(a)0,^) D , 

for <f>,if) G L 2 (Xr). Again, this definition does not depend on the choice of a. 
We shall denote by II the periodization operator 

Ua(z, b,r) = S (7 ■ z, 7 • b,r). 
The class of symbols in Theorem H] is given in: 

Definition 4.1. We denote by ITS^ the image of under II, where is defined in 

AM}. 

For a G ILS^, we can use Helgason's integral representation theorem 13.21 to have an 
alternative formula for the action of a pseudodifferential operator on the quotient: 

(4.6) Op r (a)^-(z) = / a(z,b,-iUj)e^ +v ^ x ' h %(db). 

J B 

This expression makes sense for all values Vj G C since a(z, b, r) has a holomorphic continu- 
ation to r G C. 

There is a standard relation between the Schwartz kernel K a (z, w) of Op(a) on D and the 
Schwartz kernel K^(z,w) of Op r (a) on Xp, namely 

Kl(z,w) = ^K a {z^w), 

where the kernel K a (z,w) on D is defined in fl4.ll) . One has formally 

/ K T a {z,w)f{w)Vo\{dw)= I K a (z,w)f(w)Vo\(dw) 

for any T-automorphic / and a. On the other hand, has a distributional eigenfunction 
expansion 

(4.7) K r a (z,w) = J2(OMa)<f>i,<f>k) JMfaW- 
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5. Intertwining the geodesic flow and the Schrodinger group on the 

universal cover 

In this section, we prove the intertwining formula f ll . 12[) on the universal cover D. This 
is done by defining analogues of Wigner and Patterson- Sullivan distributions on D and by 
finding an explicit relation between both. 

5.1. "Wigner distributions" on D. 

Definition 5.1. For b,b' G B and u, v' G zR, the Wigner distributions W( Uj b),(v',b>) £ V(&Dx 
R) are defined formally by: 

/ a(z, b, r)W^ b ) t ^r }b ^(dz, db, dr) = (Op(a)e U}b , e v , iV ) 

for a having the symmetry ( 14. 2\i . 

For b G B, we denote 5b{db) the distribution density on B corresponding to the Dirac mass 
at b, defined by J B f(b)Sb(db) = f{b) for every smooth /. 



Proposition 5.1. We have: W(y,b),(y',b>)(dz,db,dr) = e v ^{z)e u i^ i {z)8b{db)5-i v {dr)Vo\{dz). 
Proof. If v — ir, it is immediate from the definitions that 



(Op(a)e V)6 , e v > >v ) = / a{z,b,r)e irfi {z)e v ^y{z)Vo\{dz). 
Jd 

□ 

We define the Wigner transform of a function a G C^°(G x R) obeying the symmetry (14.21) 

by 

W : C7 C °°(C7 x R) L 2 (B x iR x B x iR, db <g> p(dr) <g> db' ® p(dr')) , 

Wa(v,b,is',b') = W^ b) ^ v ^b'){a). 

Note the "minus" sign in front of 2/. The following proposition proves completeness of the 
Wigner distributions. 

Proposition 5.2. The Wigner transform extends to L\ V {G x~R,,dg x dp(r)) as an isometry 
and satisfies the inversion formula, 

a(z,b,r) = e -^ +ir) ^ h) - [ [ e^- tr ' ){z ' b,) Wa(ir,b,ir',b')db'd P (r'). 

Proof. For r, r' G R, the Wigner transform is given by 

Wa(ir,b,ir',b') = [ a(z,b,r)e^ +ir) ^ b) e^ +ir ' ){z ^ Vo\(dz), 
Jd 

it is the Fourier transform of a(z, b, r)e^ +ir ^ z ' b * > with respect to z, evalutated at (b',—r'). 
The inversion formula and the isometry 

(5-1) \\ a \\L 2 w (Gx-R,4gxdp{r)) = 1 1 Wa(ir, b, ir , b ) \ \L 2 (BxiR + xBxi'R + ,db®p(dr)®db'®p(dr')) 
(5.2) = -| I Wa(ir, b, ir' , b')\ \L 2 (BxinxBxiR,db®p(dr)®db'®p{dr'))- 

follow from the Plancherel and inversion formulae for J 7 . □ 
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5.2. Patterson-Sullivan distributions on D. 

Definition 5.2. For u, v' G iH, the Patterson- Sullivan distribution PS( Ut b) t (- u ',b') = PS e ^ b) , e( 
associated to the two eigenfunctions e(y,b){z) = e^ +l/ )^' 6 ) and e(_ l / ) y)(z) = e^ - "')^. 6 '} is the 
distribution on ST) = x R defined by 

(5.3) PS e( M e , , „ (d&, dr) = Wb } Sb ' {db 2 e^dr. 

We use the coordinates defined in §2.21 We note that PS e ^ b)>e , b , is undefined if b = b'. 

We chose this somewhat awkward notation (writing v' instead of — v') so that the definition 
can be straightforwardly extended to v G C when we go to a compact quotient. For higher 
rank symmetric spaces, this formula was generalized by Schroder [SchDissj . He pointed out 
the fact that v' has to be replaced by —w.v' (where w is the longest element of the Weyl 
group) if one wants the "diagonal" [y — i/ G ia*) Patterson-Sullivan distributions to be 
v4-invariant. 

We now prove an analogue of Proposition 11.11 on the universal cover. Recall that S*D is 
naturally endowed with the density e^'^ Vo\(dz)db, corresponding to the Liouville measure 
on the unit tangent bundle, or to the Haar measure in the group theoretic picture ST) = G. 
In what follows we have to distinguish between distributions and distribution densities on a 
manifold, see |Ho II Ch. VI]. The choice of a preferred density allows to identify both. On the 
boundary B (endowed with the density db) we will denote S bo (b) the distribution defined by 
the Dirac mass at a point b a , and 5t, (db) = 5f, (b)db the corresponding distribution density, 
defining the linear form / i— > f(b ) on C°°(B). We recall that distributions can be multiplied 
under certain assumptions on their wavefront sets [Ho Ij Thm 8.2.10. 

Proposition 5.3. On ST), define the distribution e u>b (z,b) = e(~2 +I/ )^> 6 )# 6 (&), corresponding 

to the distribution density e v>b (z, b)e {z ~ b) Vo\(dz)db = e^ +u){z ' b) 5 b (b)dbVo\(dz). 
We have 

{e^.LlZ^) (z, 6)e<*' 5 > Vol(dz)d~b = 2ti.2^ PS e{v b) ^ i y) (dz, db) 
where i denotes time-reversal. The product on the left-hand side is well-defined for b ^ V . 
Proof. Writing e Vib in (b, b', t) coordinates, we have 

(5.4) e„, 6 (&, b', t) = e^ +u){9{l ' l ' )at -°' l) 5b(b). 
Its time reversal is thus 

(5.5) Le v , b (b,b',t) = e (-k+^^-r^) 5b (9). 
By the identity of Lemma 12.11 of §2.21 we have 

(5.6) {g (h' 1 ~b)a tl b}=t-\og 1 -^^. 
Multiplying the two distributions gives 

(5.7) (e^.teZ^) (b,b',t) = 2- 1 +»~ 17 e '*T ^ 8 b (b)5 b ,{b>). 

\b — u \- l + u - u 
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Multiplying by e {z > b) Vo\(dz)db = 4vr jf^dt, we find 

db<5t> db' -; p^ 17 ) 

(e^.teZ^) (b,b',t)4ir - ' eft = 2n.2^'- — = 8 b {b)5y{b')dbdb' dt. 

\b — o'| 2 \b — b'\ l+u - u 

□ 

Remark 5.1. As in |AZ] . we observe that Patterson-Sullivan distributions are eigendistribu- 
tions of the geodesic flow. For a distribution density on ST) we define as the pushforward 

by g l . For g = (z,b) G G, we have e u ^{ga t ) = e^^^e^^g). This implies that 

g\{e v>b {z, b)e {z >~ b) Vo\(dz)db) = eG-"*e Vtb (z, b)e {z '~ b) Vol{dz)db. 

Proposition 5.4. For u, v' G iH, we have: 

Jpq — P -t(v+V) p q 

The proof is immediate. If we extend PS e{vb)>e ,, b , to ST) x R by taking PS e ^ b)te ,, bl) ® 
8 v _-j, and if we extend the geodesic flow to ST) x R by letting 

G t (z } b,r) = (g rt (z } b) } r) } 
we have _ 2 

G\ (PS^ ® 6^) = e i ' J ^ LJi PS W{ _ v , y) ® K^z- 
In other words, for v = ir, v' = ir', 



G\ [PS e(irtb)>e{ _ irl b/) ® 5^ j = e ' » PS^,,^,, ® 

When working on compact quotients we will have to worry about a possible extension of 
these formulae to the case of complex r, and this is why we pay attention to write formulae 
that can be adapted in a straightforward manner to r G C. 

5.3. Radon-Fourier transform along geodesies. PS-distributions are closely connected 
to the Radon transform along geodesies. As reviewed in §2.2[ the unit tangent bundle ST) 
can be identified with x R: the set represents the set of oriented geodesies, and 
R gives the time parameter along geodesies. We denote by 7^ the oriented geodesic with 
endpoints b' , b. 

Definition 5.3. The geodesic Fourier-Radon transform is defined by 

n:C c (ST))^C c (B^ xK), by Hf{bf,b,r)= [ f(g(b',b)a t ) e - trt dt. 

Jr 

It is clear that TZ intertwines composition with g l and multiplication by e irt , i.e. 

(5.8) n(fog t )(b>,b,r) = e irt 1lf(b',b,r). 
By the Fourier inversion formula, 

(5.9) f(g(b', b)a t ) = ±- f Kf(b', b, r)e irt dr. 

We call "Patterson-Sullivan transform" the pairing of the family of PS-distributions with 
a test function. 



GEODESIC FLOW AND SCHRODINGER FLOW ON HYPERBOLIC SURFACES 21 



Definition 5.4. The PS-transform is defined as follows: 
PS : C™{G x R) ->■ C°°(BW xiRx iR) on G by 



PSa(z/, b, u', b') = PS^U-vip) (a^-v) = r. 77777 - / a ( #(&', b)a T , - — ^ ) e {iy+i/)r dr 

V -sr/ p — O' Jr. V A / 



The PS'-transform is related to the Fourier-Radon transform as follows: 

(5.10) PSa(v,b,v',b') = - I = Ka v _ v {b'Ai{v + 17)). 

o — o - Si- 
Using the inversion formula for TZ, one gets the inversion formula for the PS'-transform : 

Lemma 5.5. The function a is determined from its PS-transform PSa(ir, b, ir', b') by 

a(b\ b, t, R) = -e 2iRt \b - b'\ 1+2iR [ PSa(ir, b, i(2R - r), b')e- 2irt dr. 
K Jr 

Using the Parseval identity for the 1-dimensional Fourier transform, as well as Lemma I37T1 
expressing the Haar measure in the (b', b, t)-coordinates, one gets 

Lemma 5.6. 

\\ a \\h(GxR,dg®d P ) = ~ f \PSa(ir,b,ir',b')\ 2 dbdb' (^-) tanh (tt^-^J drdr'. 

71 Jb,b'£B,r,r'&H \ * J \ A / 

This formula can be compared to that obtained for the Wigner transform (15. ip . Note, 
however, that the drdr'- density is different in the two formulae. We also stress the fact that 
we do not ask a to have the symmetry (14. 2 p here. 

5.4. Operator sending the Patterson-Sullivan distributions to the Wigner distri- 
butions. If a is a function on SD ~ G, and v G C, we define the function L u a on G 

by 

L u a{g)= i a(gn u )(l + u 2 )~ { ^ +v) du. 
Jr 

In this section, we prove the following : 

Proposition 5.7. Let a G C?{G), v,v' G iR and (b',b) G P (2) . Then L_ v {a) G C°°(SD). 
Although L -^(a) is not compactly supported, the pairing PS( Vi b)(- v ',b') {L_-^(a)) is well de- 
fined, and we have 

PS {vM -s >v) {L_ 7 (a)) = 2-^-^W (vm ^, M) {a). 

A proof of this property by direct computation was given in [AZJ, in the "diagonal" case 
v = v' G iR. The proof given there could be transposed to the general "off-diagonal" case. 
Here we give an alternative presentation based on Proposition 15.31 and on the invariance 
properties of the distributions 

To check the formula, we use that a t n u = n ue ta t and the KAN (Iwasawa) decomposition, 

(5.11) n u = fc„a_i og (i +n 2)n/-( M ), 
as in fl2U). 

We now consider the action of L_-^ on the product (e^.tET/Ty) (we note that the Poisson 
density e^ Vol(dz)db is n u -invariant, so can be taken out of the integral defining L-p). 
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Proposition 5.8. (1) e u j, is n u -invariant; 

(2) L v (te Utb )(z,b) = irer U) fy(z) (it does not depend on b, in other words, it is a right-K- 
invariant function); 

(3) we have L_-^ (e u ,b-^-u',b') = ^^,b^{-v',V) ■ 

Proof. (1) is obvious, and (3) follows from (2), so we just have to prove (2). 



(5.12) £„K*)G7) = J R (l + u 2 )-^+»he r>b {gn u )du. 

We rewrite n u using (12.71) . Then we: 

• remove the right N factor since Le u>b is right- N- invariant; 

• replace the A factor inside ie Vjb by a factor e ( ~2+ i/ ) lo g( 1 + M2 ) outside, since ie Vjb is an 
A-eigendistribution of eigenvalue = | — v which is evaluated for a_i og (i+ u 2). 

• then change variables to K with 9 = arctanw. 

We then have, 



(5.13) 



L u (te Utb )(g) = J r (1 + m 2 ) l (ie U)b )(g)(gk u )du = ix j K (it Vyb )(gk)dk 



7T j K e^ b {gk)dk = ttJ e„ jb (z, V)e^ v >dbf = 7red +i/ )<*> fe > 

for g — (z, b). □ 

Multiplying by the Poisson density, we finally get 

(e U)b .Lt_ u , e {z ' b) Vol(dz)db = ne^^e^^z) Vol(dz)5 b (b)db. 

We recognize from Proposition 15.11 the expression of the Wigner distribution, for a test 
function a that does not depend on the r-parameter. Comparing with Proposition I5.3[ we 
obtain Proposition 15.71 

Proposition 15.71 was proven for a function a defined on <SD, in other words a function on 
ST3 x R that does not depend on the last variable. In the sequel, we will apply Proposition 
15.71 to an arbitrary function on SD x R, using it in the following form. If a is a function 
on SD x R ~ G x R, and r 6 R, we define the function a r on 5*D ~ G by a r (g) = a(g, r). 
Proposition 15.71 implies that 

2- (1+ ^ 7) WW) ( - i /,6 (a) = PSfuM-^v) , 
for v = ir G iH, and u' G iR. 

5.5. The operator C Recall that we have extended the PS-distribution -?>%,.,&),(-«•',&') 
(r, r' G R), originally defined on SD, to SD x R, by tensoring it by 8r±r± on the R- variable. 

2 

We now look for an operator C that acts on functions (distributions) defined on SD x R, 
with the property that PSCa(ir, b, ir', b') = Wa(ir, b, ir', b') (r, r' G R). This means that we 
must have PSCa{ir,b,ir' ,b') = 2 1+%r+ir PSW^j-ir'^') (L ir /(a r )). 

By the PS-inversion formula (Lemma 15.51) . we have for all (b', b) G B>( 2 \ t G R, R G R, 

(5.14) Ca{b',b,t,R) = - e 2iRt \b-b'\ 1+2lR / PS{Ca){ir,b,i{2R - r),b')e~ 2irt dr 

^21+^ r ^ +u 2y { l+iR) ar0h u^ }bjT yi{R-r)(t-^ 

* Jr 
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In other words, letting g = (&', b, t), 

2l+2iR r 

(5.15) Ca(g,R) = / (1 + u 2 )~^ +iR) a(ga i og(1+u2) n u , r)e 2i{r - R)l drdudr. 

7T J T 2 

5.6. Intertwining. In this section, we prove that the operator C intertwines V t and G t on 
D. We recall the Hilbert space in Proposition 15.21 

Theorem 7. The intertwining operator £ is an isometry from L\ V {G x H,dg x dp(r)) to 
the space T-Lps(D) of functions such that 



1 

4 

and we have 



\PSa(ir,b,ir' ,b')\ 2 dbdb'p(dr)p(dr') < +oo. 



(5.16) £oF* = G*o£, 

where both sides are bounded operators from L\ V {G x H,,dg x dp(r)) to "Hps(D). 



Proof. First, we consider the action of both sides on a 6 S™. We know that the space 

is preserved by V t . We shall first check that Ca is a continuous function when a £ «S^, and 

that (15 . 161) then holds as a pointwise equality between £(V'a) and G*(£a). 

For g = 6), we see that 
(5.17) 

Ca(z,b,R) = ^ / (l + « 2 )- ( ^ +i;i) « ° ^" ° r^f^^le'^rfrWr 



= t^l J(l + M 2)-(^-Hfl)(/ - a?)* a o fr« o ^- l£S ^ (z, 6, r )e 2l ^-^^r^ ITT 5 3F . 

where the r-integration by parts is used to gain powers of r and make the r-integral 
convergent. We also know from the definition of that |(J — d 2 ) N a(z, b,r)\ < Cn,m,xo{^- + 
r .2^-M e -Md(z,x ) £ Qr an y ]\[^ M > and any given x D . 

If z stays in a fixed compact set, denoting (z, b) = h u o g T ~~ j ^r -(z, b), one can check by 
hand that exp d(z, x Q ) > Ci(l+ |u|)e' T ' — C 2 , with C\, C 2 > 0. For instance for (z, b) = e £ G, 
we compute explicitly 

e T / 2 _ue^/2 
(l +M 2)l/4 (l+«2)l/4 

e-^^l + u 2 ) 1 /^ 



Q log(l + u 2 ) 
T 2 



In the Poincare upper plane model, identified with PSL(2, H)/K, this element represents a 
unit tangent vector based at 



ue 



(l + U 2)l/2 (1+ M 2 )V2- 

The hyperbolic distance of this point to the origin x Q = i is given by 



cosh d(z, i) = 1 + 2 



-u' + uVT 



4(l + n 2 ) 1 /V T >8(l + M 2 ) 1 / 2 el T l-7. 
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It follows that 

(5.18) |(/ - d 2 ) N a(z, b,r)\< C(l + r 2 )^(l + u 2 )- K/2 e~ N ^ 

(with C uniform as z stays in a compact set), so that the integral (I5.17P does make sense, 
and defines a continuous function of the variables (z,b,R). To keep this paper reasonably 
short, we do not investigate the additional regularity properties of £a. 

Using (15.181) . we also see that, for fixed (z, b, R), we have \{£a)og t ( y z, b, R)\ < C z ^^M e ~ M ^ 
(for M > arbitrary), and it follows that PS£a(ir, b, ir', b') is perfectly well defined for any 
r,r' G C, (b',b) G B^. The Wigner transform Wa(ir,b,ir' ,b') is also perfectly well defined, 
and £ has been constructed so that PS£a{ir, b, ir', b') = Wa(ir, b, ir', b'). We see that 

PSiG'Ca)^, b, ir', b') = e^^^ 21 PS{£a){ir, b, ir' , b') 

= e ~ i (r2 ~i 2)t Wa(ir, b, ir' , b') = W(y*a)(zr, b, ir' , b') = PS(£U'a)(ir, b, ir' , b'), 

and inverting this formula we get that G 1 £a = £V l a, for all a G (and this equality holds 
pointwise). 

We can now easily extend the intertwining formula to a G L 2 V {G x R, dg x dp(r)). Using 
formula (15. ip . we see (in a tautological way) that C is an isometry from 7^(G xR,rfgx dp(r)) 
to "Hps(D), and we have G l o £ = £ o V 1 , where both sides are bounded operators from 
L 2 W (G xR,dgx dp(r)) to H PS (T>). 

□ 

Comparing with Lemma [5.6[ we note that the norm on Hps(D) is not equivalent to the 
norm on L 2 V (G x R, dg x dp(r)) (the (irdr'-densities differ by an unbounded factor). 

In the Section HI we will mimick this construction to build two Hilbert spaces T^ty(Xp) 
and %ps(X r ) formed of T-invariant symbols, such that £ sends "H^(X r ) isometrically 
to 77ps(Xr), and such that the intertwining formula G l o £ = £ o V f holds between 
these two spaces. On the quotient, 77iy(X r ) will be naturally identified with the space 
of Hilbert-Schmidt operators via the quantization procedure Op r , but will not be equivalent 
to L 2 ((T\G) x R). 

Remark 5.2. Since G l preserves the variable r, Proposition 15.161 still holds if we modify 
the definition of £a(g, R) by a constant depending only on R. Thus, we have the choice of 
a normalization factor for £. We note that G f l = 1 and V*! = 1, so that it is quite natural 
to renormalize £ to have formally £1 = 1. This means dividing £a(g, R) by 

J {I + u 2 )-^ +i ^e- 2i ^ T drdudr = n J (1 + u 2 )-^ +iR Uu. 

The function ^o{s) = j_ °°(1 + u 2 )~ s du (5te(s) > |) extends meromorphically to the whole 

complex plane by /io(s) = -^w^ ^ (see p. 65-66 in |Hej ). 

The renormalized £ now satisfies PS^ rt b)^ ir >^'){£a) = l^(ir,fc), (-«•', b')( a ) if we define the 
normalized PS'-distributions by 

(5.19) PS(ir£)^- ir >V) = 7l7i I - + i I PS( irj b),(-ir>,b>)- 
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Remark 5.3. When working on the quotient, we will need the following properties of C, 

which follow from its explicit expression (I5.15p . 

First note that, if a G S^, then Ca(g, R) has a holomorphic extension to R G C 
Assume that a(z, b, r) G satisfies in addition, for every e > 0, p, all q > 0, every integer 

a, and every if-left-invariant differential operator D acting on B, a bound of the form 

d a 

-Da(m,b,r) 



(5.20) supe" |r| 

(r,b) 



dr c 



< +oo, 

Cp(G/K) 



in {13771(7") | < |}. In other words, we strenghthen the definition of by asking that a 
decay superexponential fast in r, instead of superpolynomially fast. We will denote by <SJ° 
the space of such symbols. 

Then, for any fixed g G G and R G C, the map t i-> Ca(ga t , R), originally defined for 
t G R, has a holomorphic extension to t G C. In particular, (G t Ca)(g, R) is well defined for 
ReC. 

5.7. Remark on the wave flow. Let us briefly discuss the case of the wave flow, 
e it-y/-A-i/4 ^ or alternatively, e rfv/ ~^). The corresponding quantum evolution is 



(5.21) /9*(Op(a)) = e- i V- A ~ 1 /4 Op(a)e i 'v / - A - 1 /4 =; p([/'a). 

The explicit expression of {/* is given in 



Since e l V A 1 / 4 e^ b ) = e itr e^ rjb ) and e l V A 1 / 4 e ( _ if . i& ) = e i * r e(_ ir . )fe ) for r > 0, we see 
that ?7* defines a unitary operator on L^{G xH,dgx dp(r)), and that U^W e(irb)je , y = 
e it(r-r') w for r.r' > 0. 

We also have g^PS e(ir b)>e{ _ ir , b , ) — 

e it(r r ')pS Hir b)!e{ _ ir , bl) (where g t is the unit-speed geodesic 

flow). 

It follows that (£ o ^^P^^^e,,^^,, = (5-* o C)^PS e{ir b)fi{ _ iri y) for r, r' > 0. So in this 
sense, £ also intertwines the wave group and the unit-speed geodesic flow. But because we 
restricted to positive values of r, the result is not apriori as strong as for the Schrodinger 
flow. 

A further defect is that e l * A / _A_1//4 does not preserve the Schwartz spaces C P (D), because 
yf. is not a holomorphic function on the complex plane. On a compact quotient, there is 

also a problem with the definition of e l *\/ _A_1 / 4 f or low eigenvalues (in particular, it is not 
unitary). Some of these problems are circumvented by using e lt ^~^. But they explain why we 
prefer in this article to work with the Schrodinger group e ltA . We discuss the intertwining 
of the wave group and geodesic flow further in [AZ2j . where the intertwining involves a 
modification of C 



5.8. Remark on semiclassical symbols and on an "exact" Egorov theorem. In the 

study of quantum chaos in the semiclassical regime, one often works with the flow (e lth ^), 
in the limit h — > 0. In this case, one considers semiclassical symbols also depending on 
h > : here, for instance, we could define them as functions a% G S<^, having an asymptotic 
expansion 

+oo 

an ~ ^2 ^ afc ' 

k=0 
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dk G S™, the expansion being valid in all the 5^-seminorms. In this setting, one works 
with the operators Op h (a) = f Op(a(z, b, Hr)). 

Introducing the operator M^a(z, b, r) = f a(z, b, hr), it is natural in this context to introduce 
the notations £ h = M~ l o Co M h , and V£ = M~ l o V th o M h . Note that M' 1 o G nt o M h = G*. 

In the semiclassical setting, the intertwining relation reads L% o V£ = G l o £ h . Using the 
stationary phase method, both sides of the equality Cn ° Vh a h = G l ° £>h a h can be expanded 
into powers of h, and on both sides the coefficient of h? is G*ao : this is an expression 
of the so-called "Egorov theorem", which says that V^aa has an expansion starting with 
G l a Q + 0(h), combined with the fact that L% — I + 0(h). 

Our intertwining relation is intended to given an "exact" form of the Egorov theorem : if 
we know that Cn is invertible, and have an explicit expression for its inverse, we can then 
define 

(5.22) 6^ h (a) = Op h (C^a), 

and this new quantization procedure will have the property that 

e- m $0{f h (a)e m $ = 0^n(G l a). 

Such an exact intertwining relation is often called the "exact" Egorov property, and is so far 
only known in the euclidean case, where the Weyl quantization Op^ 7 has the property that 
e - tth — Opjf ' (a)e lth ~ = Op h (G t a) (where a is a function on T*H d = H d x R d with reasonable 
smoothness and decay properties, A is the euclidean laplacian on R d , and G l is the euclidean 
geodesic flow). 

Being able to compute £ _1 amounts to computing C~ l , and this can be done formally 
as follows : we must have W(i r> b') j (— ir',b') (C- X a) = PS 

(ir,b),(—ir' ,b')( a ) > an d we can recover the 
expression of C~ l a using the inversion formula Proposition 15.21 We find 

L- l a(z, b, r) = e -(£+*0<*,6> [ ^ > f'Tj e^ ir >e^- ir '^ b 'Ub'dp(r')dr. 

v ' ' > I \h — h> \l+tr+ir' ' 

Jb'£B,r'>0,T£-R \° ~ I 

To express this in group theoretic terms, let us consider the special case z — 0, b — 1 (in the 
disc model), corresponding to g = e G G. Using the calculations of §2.4[ we find 

// r _|_ r '\ -l + ir + ir' . 2-( 1 + ir+ir ') 
a I n u a log (i +u 2) ; ] (l + u 2 ) 2 e ( ir ~ ir > T dudp(r')dr. 

\ T + 2 2 J 7T 

More generally, using G-equivariance of the formulae, 

// r _)_ r '\ -i +ir+ir i 2-( 1+ir+ir ') 
a \gn u a_ | log(1+u 2 ) ; — - — J (l + u 2 ) 2 e i"— ir n _ dudp(r')dr. 

Using f)5.22p . one defines a new quantization procedure satisfying the exact Egorov prop- 
erty. A drawback is that the regularity properties of C^ 1 are not well understood. One can 
also note that, contrary to what is usually expected from a quantization procedure, Op^(a) 
is not the multiplication by a if a = a(z) is a function on G/K. 
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6. Intertwining the geodesic flow and the Schrodinger group on a 

compact quotient 

6.1. A remark on the spectrum of the geodesic flow and the Schrodinger flow. In 

this section, we investigate the meaning of the formula 

(6.1) CoV^G'oC 

on a compact quotient. We already saw that this formula holds as a pointwise equality, when 
a G S™. If a formula such as (16 .ip holds on a compact quotient, one would be tempted to 
infer that there is an explicit relation between the spectrum of the geodesic flow and the 
spectrum of the Schrodinger flow. On the other hand, one might object that the spectrum of 
g l on L 2 (r\c7) is continuous (a manifestation of the ergodicity of the geodesic flow), whereas 
the spectrum of the quantum evolution a 1 on the space of Hilbert-Schmidt operators is 
discrete. Of course, there is no contradiction at the end, since the spectrum depends on the 
space where an operator acts, and since £, £ _1 are unbounded operator. 

In fact, the theory of resonances for the geodesic flow already describes its spectrum on 
spaces other than L 2 , and finds some Banach spaces of distributions where the geodesic 
flows has some discrete spectrum (see the work of Ruelle and others, e.g. |BT[ IBKL] |BL[ 
IFRSl RjLl ILivt lRu87l Rugh92 Rugh96|). This approach was especially developed to describe 



the correlation spectrum of the geodesic flow for smooth (or Holder) functions. As recalled 
in [AZj . the distributions e u . are the generalized eigenfunctions arising in the resonance 
expansion of the geodesic flow. We must also note that the e v . are precisely the off-diagonal 
Patterson-Sullivan distributions PS U .^, associated with the pairs of eigenfunctions 4>j an d 

the constant function = 1, with spectral parameters, respectively, Vj and |. As to the 
other off-diagonal Patterson-Sullivan distributions PS Uj - Uk , we saw they are also generalized 
eigenfunctions of the geodesic flow, but their spectral interpretation remains unclear : they 
do not appear in the usual resonance theory. 

In the next section, we shall construct a Hilbert space of distributions Hps, on which 
the geodesic flow acts, with dual eigenbasis the whole family PS v - Vk . In Definition [TJ we 
introduced the Patterson-Sullivan distributions PS Vj _ Vk = PSy,^),^-^) associated with 
the pair of eigenfunctions (0^,0^) and the choice of spectral parameters (uj, —i>k)- On the 
quotient, the space Hps will be constructed by transposing the analysis done in §5J replacing 
the distributions PS e{vb)te , b , by the family PS Vj _ Uk . A drawback is that elements of Hps 
cannot be characterized in a simple way in terms of their regularity/growth properties (for 
instance, Hps does not coincide with one of the known Sobolev spaces). However, we will 
give sufficient conditions for a distribution to belong to Hps in Section [71 

Proposition 6.1. The distribution PS Uj - Uk is T-invariant. 

Proof. Recall that 



PS^ k (db',db,dr) = l_ n ;tJ ^ + ^dr. 
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We use the formulae of Section [3j and the fact that 7^7^ (do) = e-^ +u ^°'^T v .(db) for 

7 g r. 

77 1 PS (db' db dr) = ^ (fgTf^gy) ( , J+PE)(r+ ^>W) )rf 

= c -(h+Vi)h-o,Tb) c -a-Vk)(TQ,TV) Tu 3 ( db ) T -"k W g l (l + ^-Pfe)(( 7 . , T b) + ( 70 ., 7 .6')) e (^+Pfe)( T+ <7-°.7-fc>-<7- ,7 

16 — i/ll+^i— "fc 



T Vj {db)T- Vk {db ) , 

______ e K u rfr = PS Uj< _ Uk (db,db,dr). 

□ 



6.2. Hilbert-Schmidt symbols. We first have to remark that it does not make sense to try 
to generalize Proposition (15.161) to T-invariant symbols, in the form of a pointwise equality 
£r ° V^a(z, b, r) = G\. o £ r a(z,b,r). Indeed, really acts on symbols via the definition 

e - lt ^r Op r (a)e u ^r = Op r (V», but on a compact quotient there is no inversion formula 
allowing to recover the symbol a in a unique way from the operator Op r (a). In fact, one 
can easily see that two different symbols on SXr x R can yield the same operator on the 
quotient. 

To put this another way, consider two symbols a, b G which have the same periodiza- 
tion a = f Ha = lib. We can try to define V^a by VpCt = UV t d, but there is no reason 
why we should have UV t d = UV t b pointwise, in other words this definition of V^a will de- 
pend on the choice of a, b. On the other hand, Op r (V^a) does not depend on any choice, 

and coincides with e ~ Op r (a)e lt ~ . This explains why we shall not try to prove that 
£r ° V^a(z, b, r) = G r o C T a(z, b, r) pointwise - one needs in some sense to quotient out by 
the symbols a such that Op r (a) = 0. 

The aim of the section is to present a (simple) construction of two Hilbert spaces Hw and 
Hps, having respectively the families (W v .- Vk ) and (PS U] ^ Uk ) as dual orthonormal bases, 
with the following properties : V 1 acts unitarily on %w, G 1 acts unitarily on Hps, C sends 
Ww isometrically to Hps, and the intertwining relation (16 .ip holds on these spaces. This 
constructions mimicks what was done at the end of §5.61 with the additional difficulty just 
mentioned, of having to quotient out by the symbols a such that Op r (a) = 0. 

We start with the Hilbert space HS(Xr) — L 2 (Xr x Xr) of Hilbert-Schmidt operators 
on the compact quotient Xr = T\G/K. On this space, the quantum evolution a 1 has the 
orthonormal spectral expansion 

(6.2) a 1 = (<f>, ® «) ® (<f>j ® <f>D* 

j,k 

(here, v\ G R, but we wrote Vk 2 to underline the fact that the dependence of the expression 
is antiholomorphic w.r.t. (pk^k)- The Hilbert-Schmidt norm is defined by ||^4||#5(x r ) = 
Tr(yL4t), associated with the scalar product (A,B)hscx.t) = Tr(ABt). Starting with a T- 
invariant symbol a belonging to n^^ 3 , we obtain a Hilbert-Schmidt operator Op r (a) G 
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HS(X-r), with norm 

II Op r (a)||L(x r ) = Tr L2(Xr) Op r (a)Op r (a) t 
= ^|TrOp r (a)0,®0*| 2 = ^|W^(a)| 2 . 

j,k j,k 

This suggests to define the Hilbert space "Hyp as as follows: 
Definition 6.1. Tiw is the completion of the symbol space HS^ with respect to the norm 

The scalar product on 7i w is defined by 

^b) w = Y^Wl k {a)Wyb). 

Here we have to stress two important facts : 

• The norm ||.||w is actually a seminorm. Again, this comes from the fact that we can 
have Op r (a) = although a ^ 0. The reader should not be too surprised by this 
fact, which already occurs in the euclidean case when one wants to study the Weyl 
quantization Op w (a) of a symbol a(x,£) (where G T*R d ), 27rZ d -periodic in the 
x- variable. If a ^ 0, then Op w (a) defines a non- vanishing operator on L 2 (R d ), but 
the periodization of this operator may vanish when acting on the torus L 2 (R d /27rZ d ). 
Actually, this happens if a(x, £) vanishes when £ is a half-integer. In the hyperbolic 
setting, a similar phenomenon occurs. Although the operator Op (a) is non-zero on 
L 2 (D), its periodization Op r (a) can vanish when acting on L 2 (Xp). The difficulty is 
that there no easy characterization of the symbols a such that Op r (a) = 0. 

• Another related issue is the following. On the universal cover D, we have seen that 
the Hilbert-Schmidt norm of Op(a) coincides with the L 2 -norm of a, seen as a function 
on ST) x R endowed with the measure e^'^dzdb x dp(r). This is no longer true on 
the quotient. In other words, for a T-invariant symbol a, \\a\\w is n °t the L 2 -norm 
of a on SXp x R (see §7.2p . Again, the same phenomenon already occurs in the 
euclidean case. 

By definition, Wj k is a bounded linear functional on Hw for all j, k, in other words 
Wj k G Hyy-i where i s ^he dual Hilbert space to T-Lw The Riesz theorem endows W^r 
with a dual inner product, and the Wj >k form an orthonormal basis of H^- For b G T-Lwi t ne 
series 

(6-3) E^WliMW^ 

converges in W^. In fact, the operator W? fc («)W£ fc : Ww T~i-w i s J us ^ the standard 
(antilinear, unitary) isomorphism b — »■ from 7i w — > 'K* w - 

As described earlier, we can define Vp acting on HS^ by V^a = IlV t d, if a = lia and 
a G S™. This definition depends on the choice of a, however W^ k {V^a) does not. We 
then note that the evolution Vp can be extended to %wi an d is obviously unitary, since we 

2 2 

have Wj k {V^a) = e ftJ T— W^ fc (a). The Wi gner distributions form an orthonormal basis of 

2 2 

eigenf unctions of the adjoint V^ in "H^ : we have V^{Wj k ) = e lt 2 Wj k . 
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It is difficult to find a full characterizations of elements in Tiw i n terms of usual Sobolev 
spaces. However, in £J7J we shall give a sufficient condition for a function to belong to of Tiw, 
in terms of its regularity and decay rate at infinity (in the r variable). 

Carrying this forward one step, we define another Hilbert space Tips- We consider the 
space n e C(R e , C£°(G)) of symbols a(g, r) which are holomorphic in r G C, taking values in 
the space of smooth compactly supported functions of (z, b) (here the extension to r G C is 
needed to deal with the case of low eigenvalues of the laplacian); and such that, for every 
integer a and any k, q > 0, for any e, we have 

Qa 

'! + \r\) Q TT-^a < +00. 



sup 

r£R E 



dr c 



c k , 

z.b 



Definition 6.2. We define Tips as the closure of II (n e C(R e , C%°(G))) under the "scalar 
product" 

(6-4) {f,g) PS = J2 PS l-^f) PS ^M- 

An alternative which is closer to Definition 16.11 would be to use the normalized Patterson- 
Sullivan distributions ( I5.19p . Here we have to choose a value of the spectral parameter Uj 
for each j, and we make the standard (but somehow arbitrary) choice : Vj G [|, 1] U iR+. 
The corresponding (semi)norm will be denoted H-Hps- We then have PS„ _ Vk G Ti PS , and 

k PS^._ Uk (»)PSl _ Uk : Tips — » Tips ls ^ ne standard (antilinear, unitary) isomorphism 
from Tips 7~P PS . The evolution Gf can be extended to a unitary operator on Tips, or to 
Tips by duality. 

6.3. Proof of Theorem [4] . We first state a more complete and precise version of the 
theorem. The Hilbert spaces Tiw, Tips are discussed in more detail in 

Theorem 8. We have: 

(1) Vp is a unitary operator on Tiw and on Ti* w . The Wigner distributions form an 
orthonormal basis of eigenvectors ofV^ in Ti* w . 

(2) The maps £p '■ Tiw — > 'Hps, an d £rtt : Ti* PS —> Ti w are isometric isomorphisms; £rj 
sends PS Uj - Uk to Wj >k - For a G HS^ , we have 

£ r o V^a = G r o C r a, 

as an equality between two elements of Tips- 

Proof. The intertwining relation on the quotient follows from the intertwining on the uni- 
versal cover, and from the relations 

Wl k {a) = W hk {a) = j ^ ^^ ^(^^.(^^(^O, 

(6.5) PSl^ Vk {a) = PS Vj ^ k {a) = J PS^,^ bl) {a)dT v .{db)W^ k {db'), 

valid for any T-periodic a, and any a such that a = Ua, provided a has the sufficient 
smoothness and decay so that these expressions are well defined. 
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Take a £ ITS^ 1 , that is, a = Ua with a £ S™. We know from the calculations in §5] that 
PS^iJCa) = J P^, 6) , e( _^ 6/) (£S)^.(rf6)rf^:(^) 



Theorem [3] is proved this way, and we also get a proof of Proposition II. 1\ by combining 
Proposition 15.31 on the universal cover, with formula (16. 5p . 

Now, take a £ LLSJ , that is, a = Ha with a £ (see Remark 15. 3p We know, from 
Remark [5. 3 1 that G*£a(g, R) is well-defined for all R £ C, and we have 



2 2 

"3 ~"fc 



e«— s-W^a) 
W j , k (y t a) = PS v ,. Vk (CoV t a) 



The identity holds for all j, k (everything is well defined even for low eigenvalues of the 
laplacian) . 

From this, we can deduce the following : 

• we can define C^a £ "Hps by PS^ _ Vk (Cra) = PS Vj - Vk {Ca), and this definition does 
depend on the choice of a. 

• the adjoint £rjt sends PS^._ Vk to Wj k , hence £rjj is an isometry from Wp S to H^, 
and £r is an isometry from Hw to Hps- 

• the family PS^._ Vk forms an independent family in H* PS : if Yljk a 3kPSy._ Uk = 
with ^2 \o>jk\ 2 < +oo, then atjk = 0. This comes from the fact that the Wj t k form 
an independent family in W^, and Cr$PS$._ Vk = Wj k . It follows that the family 
PS^_ _ Uk is an orthonormal basis of Wp S . 

• for a £ ILS^ , we have 

£ r o V£a = C7p o C r a, 
as an equality between two elements of Hps- 

□ 

We note that Hilbert-Schmidt symbols are a very special class of symbols. But it is mainly 
a condition on the r-growth of symbols. By multiplying by powers of r, the intertwining 
operator and the time evolution extend readily to more general symbols. 



7. Further discussion about the Hilbert spaces H w and Hps 

We now discuss the elements of the Hilbert spaces Hw and Hps in more detail : we 
describe sufficient conditions for a function to belong to Hw and Hps, in terms of regularity 
and decay. We use the regularity properties of the boundary values T v . of eigenfunctions, 
described by Otal [D]. In the automorphic case, the regularity properties can be read off 
directly from the automorphy equation (see e.g. [MS] IMS2] ). 
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7.1. Holder continuity of T v . Following Otal [0], we say that a function F denned on R 
is 27r-periodic if there is a constant C such that F(x + 2ir) = F(x) + C, for all x. If F is 
locally integrable, its derivative DF yields a well-defined distribution on S 1 = R/27rZ, 

DF(cp) = - jf ^ ^F(0)d0 + ip(Q) [F(2n) - . 

for every smooth function if on S 1 . For < 5 < 1 we say that a 27r-periodic function 
F : R ->■ C is <5-H61der if \F(9)-F(9')\ < C\9-9'\ s . The smallest constant is denoted \ \F\\ S . 
We denote the Banach space of 5-Holder functions with norm \ \F\\s by Ag. 
We recall: 

Theorem 7.1. f[0] ) Suppose that (f) is a laplacian eigenfunction of eigenvalue — s(l — s) = 
— (| + r 2 ) , twf/j s = \ + ir and 9fte(s) > 0. Assume that \\<j)\\oa < 00 1 1 V0| |oo < oo. 
T/ien zfo Helgason boundary value T s ^ is the derivative of a Ke(s) -Holder periodic function 
F. 

In addition, letting 5 = 3f£e(s) ; we have 

\\F\\ S < t^tM (Halloo + IIV0IU), 

where C > is an absolute constant, and C s = J + °° Jq W e~^ 1+s -'*P s (tanh |, 9)d9dt; where P 
is the Poisson kernel of the unit disc. 

Outside of the finite number of "small eigenvalues" of Xp, we have 3?e(s) = | and hence 
Ty. is the derivative of a Holder 1/2-continuous function. The upper bound on ||.F||,5 given 
by Otal's proof is quite crude, but will be sufficient for our purposes. 

The behavior of C s for s = | + ir and r — > ±oo can be evaluated by the stationary phase 
method. The calculation is routine but not completely straightforward because the domain 
of integration is non- compact. However, for the sake of brevity we omit the details. One 
finds that C s ~ Cr~ 1 / 2 , with C ^ 0. 

7.2. Hw an d L^(GxR, dgxdp(r)). In this section, we clarify the relation between Hilbert- 
Schmidt inner product, which induces the inner product (,)w on symbols, and Ly V (G x 
R, dg x dp(r)). The second term in the following proposition is the discrepancy between the 
|| • \ \w and the L 2 norm on symbols (again, we stress the fact that this discrepancy would 
also appear in a euclidean situation). We denote T> a fundamental domain for the action of 
r on D. 

Proposition 7.2. Let Op r (a) be a Hilbert- Schmidt pseudo-differential operator on Xr with 
complete symbol a. Then IHI 2 ^ = f ||Op r (a)||| f5 / Xr ) is given by 

\\a\\ 2 w = j v j B j R+ \a{z,b,r)\ 2 e^Vol{dz)dp{r)db 

+ E 7e r\ W I eP , (MeBx R + fo > r)a( 7 ■ ^^e^^e^^dp{r)db V6L(dz). 

Proof. Recall that K^(z, w) = K a (z, ■yw), and that K a is invariant by the diagonal action 
of T : K a (yz, 7w) = K a (z, w). 
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The composition formula for the kernels is : 

K$oKZ(z,w)= [ K a (z,v)Kl(v,w)Vo\(dv). 
Hence K v a o K^(z, w) = J D K a (z, v)K£(w, v) Vol(dv). Taking the trace. 

Tr(K T a oK T b ] ) = Y, I I K a {z,v)K h {7 ' z,v)Vol(dv)Vol(dz). 

7gr JzeV JveD 

The rest of the calculation proceeds as in Proposition 14. 1\ using the Fourier inversion formula. 
The first term corresponds to 7 = e and the second term to 7 ^ e. 

□ 

7.3. Hw We now describe a large class of elements of the Hilbert space %w That is, we 
determine sufficient conditions on a so that Op r (a) is Hilbert-Schmidt, or in terms of Wigner 
distributions, so that 

(7.1) ||a|| 2 w = ^|^(a)| 2 <oc, 

In the following, (x) = (1 + Ix) 2 ) 1 / 2 . If C is an operator, we define ad(A)C = [A,C]. We 
denote Xj = — (f + rf) the laplacian eigenvalues. The following is of course not optimal, 
but gives an adequate idea of a large class of elements in T-L w . 

Proposition 7.3. If sup rgR+ (r) 6 | \a r \ \c& < 00, then a G 7iw 

Proof. It suffices to prove the following 

Lemma 7.4. Let 

1 1 |o| 1 1 = sup sup (Xj) 2 \\ (/ - Y 2 ) A 2 z a\ + Tj \ (I - Y 2 )A z V z a\ + r 2 | (/ - Y 2 ) V 2 a| + r 2 | (/ - Y 2 ) V z 
If 1 1 \a\ 1 1 < oo ; then a G 7iw- 

To prove this, we first note that, by Weyl's law, J2j k(^j)~ 2 (^j ~ ^k)~ 2 < 00 in dimension 
two. 

We will also use the expansion of e v . into i^-Fourier series, which takes the form 

meZ 

with Y(pj tm = 2im(j)j trn . We use the fact that ||0j,m|U 2 (r\G) = 1, proved in [Z] (the full 
definition of (fij im can be found in Proposition 2.2 of [Z], in particular, JiO = 4>f). 
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Let us write B = Op(6) = (ad(A) 2 Op(a)) o A 2 . Then, 

= E mez (2m)- 2 (^,(/-n^, m ) 

= E mG z(2^)- 2 ((/-n^- ; 0A,-) 

< sup |(/ - Y 2 )b rj \ E mez (2^)- 2 (I0,I, |&,* 



< CsupKz-y 2 )^.,, 

where C is a uniform constant. Here we use that (\<f>j\, \4>k,m\) < 1 by the Schwartz inequality 
and the fact that ||0j||l2 = ||0fc,m||L 2 = 1- 
It follows that 

\WUb)\ < C (A,)~ 2 (A, - A fc )- 2 su P |(7 - Y 2 )b rj \ 

<C(A,)- 2 (A j -A fe )- 2 

sup (A,) 2 [|(7 - Y 2 )A 2 z a\ + r,|(7 - Y 2 )A z V z a\ 

(z,b)ev 

+ r 2 \(I-Y 2 )V 2 z a\+r 2 \(I-Y 2 )V z a\] 

Here, we use that the complete symbol of Op(a) o A 2 is (| + r 2 ) 2 a(z,b,r). Further the 
complete symbol of ad(A) Op(a) is given by A z a + (| + ir) V^a.V z (z, 6). □ 

7.4. The Hilbert space %ps. We now consider the analogous question of conditions on a 
so that 

(7.2) IHI^^El^,-^ )!^ 00 - 

Proposition 7.5. 7/ a zs T automorphic and sup r (r) 12 | \a r \ |c 3 < +oo, £/ien a e "Hps'. 

This follows from 
Lemma 7.6. PFe /iov e, /or any M, 

||a||^s < Yl l^| 3/ >fe| 3/2 (ll^lloo + l|V0,||oo)(||0 fc ||oo + IIV^IU)^ +W)~ M 

3,k 

sup \b - b'\- (1+ ^- w y)n(dt) M a,_ 1 ^ K {b', b, i{vj + Pfc)), 

(b',b)&BxB L 2i 

9 -|6 - Z/l-^-^)^)^^^, 6, <(!/,- + 1*)), 



96' 
96 



^|6-6r(^-^))7e(^) M a^(6',6,i(^+^)), 



, ' 2 7 | 6 _ 6 '|-(l+^- fc ))^(^)M a ^ (6 / )6) , (z/j+ ^ )) 



9696' 
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Proof. We use the relation 

(7.3) PS T ^_ Vk {a) = J - _ * j - - Kxa*j£g$, b, ifc + ^))T v .(db)TL7 k (db'), 

which is obtained from f)5.10p . Since X a^-i^ is compactly supported on G, then the Radon- 

2i 

Fourier transform IZa is compactly supported in the variables (&', b) G B^, so the singularity 
of 1 6 — b'\ on the diagonal is not a problem. It follows by repeated integration by parts in dt 
that if a e C C M (C7), then Ha(b, b', ifa + = 0((j/j + v£)- M ). 

Let us call F v . the Holder function such that T v . = F' , in the sense of §7.11 We use the 
formula 

j tp(b',b)T Vj (db)T_ Vk (db>) = <p(0, 0) [F Vj (2tt) - F Vj (0)] [F_ Vk (2tt) - F_ Vk (0)] 

-[F Vj (2*) - F Vj (0)] f ^(b',0)FZ;(b')db> 
~[F_ Uk (2n) - F_ Uk (0)] j ^(0,b)F Uj (b)db 

+ / Jj^rtb\b)F Ui (b)FZ- k (b')dbdb', 

valid for every smooth function (f on B x B. 
It follows that 



f(t' 1 t)T,,WT-„«| 

<l|FJfe||F^| k ^ B (|^,6)|J^,6)|J^,»)|,|^(V,6)|) 

where 8j — \ + dte(vj), and the Holder norm is the one appearing in Theorem 17.11 
We can then write 

\PS^-„ k (a)\ < (v j +Vk)- M \\F Vj \\ ii \\F- Vk \\ Sk 



sup 

(fe',6)esxs 



| 6 _ br il+Uj - Vk) )n(dt) M xa^(b', b, %{v 3 + ^)), 



|_| & _ b'\- {1+v '- Vk) )n(dt) M X a^{b\ b, i( Vj + 

OO 2i 

8 \ b _ b'l-^-^Mdt^xa^ib', b, i(uj + v$), 



db 



<P -\b - br {1+u *- Uk) )K(dt) M xa^(b', Mfa + 



dbdb 
By Theorem 17. 1\ 

IK-lk = oCh-Mdl^-lloo + llv^iu 
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□ 



Moreover, by the well-known local Weyl law estimates, ||0j||oo = 0{\uj\^) and ||V0j| 
We find 

\PS Vj ,-v h (a)\ < (uj +vj:)~ M \v j \ 3 \v k \ 3 m.ax(\v j \, \v k \) \\a ^-ug \\ C M+i. 

2i 

Using the Weyl law in dimension 2, \Xj\ ~ Cj, one sees that the series 

/ , \ -N 

(rj-rfc) [rj-| |r fc | max fl^l, |r fc |) ( — 
converges for M > 1 and iV > 11. The result follows. 

We stress again the fact that there is nothing optimal in this upper bound. 

8. Appendix 

In this section we sketch the proof of Theorem [5j We closely follow the proof in Section 4 

of HZ]. 

By the generalized Poisson formula and the definition of Op(a), 
(8.1) 

(Op r (a)0 ir ,, (j) irk ) = [ ( [ X a(z, 6) e (l+^)^ b ) e (|+^)(^ b ') Vo\(dz)) T irj (db)T_ n . k (db>). 
Jbxb \Jd J 

Here we are only interested in real values of rj,r k , since we consider the asymptotics 

Tj — > +oo and \rj — r k \ bounded. We apply stationary phase to the simplify the inner D 

integral. More precisely, in |AZ] and in this article, we rewrite the integral in the form 



(Op r (a)0 ir .,0 ir J = 2( 1+ ^ +<r *> J L irkX a(b',b,T)PS ir]t . trk (db',db,dr), 

as was shown in Theorem [3l and we then replace L irk xa{b', b, r) by its expansion into powers 
of r^ 1 , obtained by the method of stationary phase. 

There is one detail that we did not discuss in |AZj . and that was mentioned to us by 
Michael Schroder (see [SchDissj ). The PS-distributions have a singularity of the form \b — 
y\-{i+irj+ir k ) on diagonal {V = b), and thus can only be integrated along functions 
that vanish on a neighbourhood of the diagonal. The function L irk xa does not satisfy this 
condition, and it is for a very special reason that its integral along PSi Tj -i rk can be defined : 
its singularity exactly cancels with \b — b'\~( 1+lVj+tVk \ However, when replacing L irk xa by its 
stationary phase expansion, one would have to justify the fact that each term, including the 
remainder term, can be integrated along PS ir . _ irk . This is not easy and wasn't discussed in 

It is actually simpler to carry out the localization step away from the diagonal with the 
original inner integral (18. ip . We see that the function 

(8.2) / X a(z,b)e^ +ir J ){z > b) e^ +irk){z > b,) Vo\(dz) 



D 



is integrated along T ir .(db)T_ irk (db'), and with the latter there is no issue on the diagonal. 

The critical set in the oscillatory integral (18. 2p occurs where V{z,b) = —V(z,b'). So 
z G yv,b- There is a neighbourhood V of the diagonal such that 7^ does not intersect the 
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support of x a f° r (b', b) G V. We take a smooth function / on B x B, supported in V, that 
is identically 1 on a neighbourhood of the diagonal, and divide the integral (18 .ip into 

/ fip', b)( [ X a(z, b)e^+ ir ^e^ +ir ^ z ^ Vo\(dz)) T ir .{db)TL~{M) 
Jbxb \Jd J 

+ [ (1 - f{b', b)) ( [ X a(z, b)e^ r ^e^ r ^ Vo\(dz)) T ir .{db)TL~{dV) . 
Jbxb \Jd J 

For the first term, the phase has no critical point, and we integrate by parts using 

Since Tj r . , T_j rfc have polynomial bounds in rj,r^, repeated partial integration shows that 
this first integral is 0((rfc) _0 °). 

The second term, because of the cut-off (1 — /(£/, 6)), is now supported away from the di- 
agonal, and can be rewritten as 2 { - l+tT ^ +tTk " 1 /(l — /(&', b))L irk xa(b' , b, T)PS irj - irk (db' , db, dr). 
The proof of Section 4 in |AZ] now applies without problem. 
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